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Summary 

An unsteady, three-dimensional Euler analysis technique 
was used to compute the flowfield of an advanced propeller 
operating at an angle of attack. The predicted blade-pressure 
waveforms for an angular inflow of 3° were compared with 
wind tunnel data at two Mach numbers, 0.5 and 0.2. For an 
inflow Mach number of 0.5, the predicted pressure response 
is in fair agreement with data: the predicted phases of the 
waveforms are in close agreement with data, whereas the 
magnitudes are underpredicted. At the low Mach number of 
0.2 (takeoff), the numerical solution shows the formation of 
a leading-edge vortex which is in qualitative agreement with 
measurements. However, the highly nonlinear pressure 
response measured on the blade suction surface is not 
captured in the present in viscid analysis. 

Introduction 

Advanced propellers are highly loaded, variable pitch 
turboprops designed to operate at high speeds and to achieve 
a higher propulsive efficiency than is achieved by the current 
high-bypass turbofans. The thin, highly swept blades of 
these propellers can produce a complex flowfield with 
leading edge vortices, tip vortices, and/or shock waves, 
depending on the operating conditions. Detailed measure- 
ments of the blade-surface steady (ref. 1) and unsteady 
(ref. 2) pressures have been carried out to understand the 
nature of the propfan flows. These tests were carried out on 
a large-scale 9-ft-diameter single-rotation, SR7L propeller 


in a transonic wind tunnel in Modane, France. The propeller 
was operated in a two-bladed configuration (fig. 1) because 
there was not enough power to run the eight-bladed design 
propeller. 

The blade surface pressure was measured on a specially 
designed, pressure-tapped SR7L blade. Unsteady blade- 
surface pressure data were obtained for propfan operation 
with a 3° angular inflow and a wake inflow, over a range of 
Mach numbers from 0.02 to 0.7 (ref. 2). These pressures 
were measured at two radial stations (fig. 2) with pressure 
transducers on both suction and pressure sides. At the 
inboard station (r/R = 0.641, where r is the radial distance 
to the point and R is the blade tip radius), there were seven 
pressure transducers on each side of the blade, while at the 
outboard station (r/R = 0.906), there were six on each side. 
Tests were conducted at static propeller conditions, take- 
off conditions, and different power-loading conditions at 
Mach 0.5 (including one at the design-cruise power loading, 
on a per-blade basis). 

Comparisons of Euler solutions and Modane wind tunnel 
test data for steady flow (0° inflow angle) conditions were 
made over a range of Mach numbers by Nallasamy et al. 
(ref. 3). The present investigation evaluates the ability of an 
unsteady Euler analysis technique to predict the unsteady 
pressure distribution on a propfan with angular inflow. It 
compares the predicted pressure-time histories for an angular 
inflow of 3° with the Modane wind-tunnel test data for two 
inflow Mach numbers, 0.5 and 0.2. 



Numerical Solution of Unsteady, Three- 
Dimensional Euler Equations 

The unsteady, three-dimensional Euler equations govern- 
ing the in viscid flow through a propeller are solved employ- 
ing a solution procedure developed by Whitfield etal. 
(refs. 4 and 5). In this procedure, the Euler equations in 
conservative, differential form are transformed from a 
Cartesian reference frame to a body-fitted, curvilinear ref- 
erence frame. Then the transformed equations are dis- 
cretized employing a finite volume technique. An 
approximate Riemann solver is used for block interface flux 
definitions, and a lower-upper (LU) implicit numerical 
scheme is used to solve the discretized equations. (Further 
details of the solution procedure may be found 
in refs. 4 and 5.) The flowfield is represented by multiblock 
composite grids to limit the core memory requirements. 
When the solution at each time step is updated, only one 
block is stored in the core memory while other blocks are 
stored on solid-state storage devices (SSD). The authors 
also used this solution procedure to compute the unsteady 
flowfield of an eight-bladed SR7L propeller with angular 
inflow at cruise conditions (Mach 0.8) (refs. 6 and 7). The 
present computations were carried out on the Cray Y-MP 
computer at NASA Ames Research Center. With the time 
step and grid resolution employed in this study, a complete 
run for one test condition took about 33 CPU hr. 

Computational Grid 

The configuration considered is that of the two-bladed, 
SR7L Modane test, which is shown in figure 1 along with the 
direction of rotation and the axes of reference. In the 
computations, azimuth angle <j> was measured from the 
z axis, as shown. This azimuth angle (<(>) reference is used 
only in figure 4; in all other figures <> =0 corresponds to top- 
dead-center, as in the presentation of experimental results in 
reference 2. An H-grid is employed to represent the flow- 
field. This grid was generated by a special purpose mesh 
generator developed for propfans (ref. 8). 

Each of the two blade passages is described by 107 by 41 
by 45 (axial by radial by circumferential) grid points, and 
each passage grid is divided into 1 1 blocks (107 by 41 by 5 


grid points in each) for computational convenience, as 
mentioned earlier. Thus 22 blocks of grid describe the 
entire flowfield (two-blade passages) with 394 830 nodal 
points. Each blade surface is represented by 49 by 27 
(chordwise by spanwise) grid points with higher resolution 
near the leading and trailing edges, the hub, and the tip. In 
figure 3, parts (a) to (c) show the bladewise surface, 
streamwise surface, and spanwise surface views of the grid, 
respectively. Figure 3(d) shows the distribution of grid 
points on the blade surface and around it. The far-field 
boundary is three blade radii from the blade tip, the inflow 
plane is two blade radii upstream of the spinner, and the exit 
plane is two blade radii behind the blade. These boundary 
locations have been found to be adequate for angular inflow 
computations (refs. 6 and 7). 

Results and Discussion 

Two representative unsteady test cases with 3° angular 
inflow were chosen for detailed comparison of the predic- 
tions and experimental data. The test operating conditions 
are shown in table I. For test 8, the advance ratio and power 
coefficient, on a per-blade basis, were nearly equal to the 
cruise values, so the nondimensional loading and flow 
angles were approximately preserved. Test 6 represents 
high-power takeoff for which the measured steady pressure 
distributions indicate the presence of a leading edge vortex 
(ref. 1). 

Hie unsteady Euler solutions were obtained for these two 
tests with the grids generated using the nominal blade setting 
angles at 75-percent radius (P 3/4 ) of the experiments (table I). 
No attempt was made to match exactly the measured and 
predicted total power coefficients. The solutions were 
carried out from an impulse start for three complete revo- 
lutions of the propeller to obtain a reasonably accurate 
solution. The results of the third revolution are stabilized as 
indicated by a periodic variation of per-blade power coef- 
ficient (fig. 4) during the second and third revolutions of the 
propeller. Running the solution through the fourth revolu- 
tion produced no recognizable change in the variables of 
interest. The results of the third revolution are analyzed, and 
pressure waveforms are compared with the experimental 
data. 


TABLE I. - OPERATING CONDITIONS FOR UNSTEADY BLADE-SURFACE 
PRESSURE MEASUREMENTS 


Modane 

Mach number. 

Advance ratio, 

Blade speed. 

Blade setting angle at 

test 

M 

J 

rpm 

75 -percent radius, 

B 3/4- 

deg 

6 

0.2 

0.881 

1671 

31.6 

8 

.5 

3.062 

1193 

55.5 
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The pressure coefficient is defined in this paper as 

P.- P 

<*-E 'tfTA 

ij O.5p_(vi + r 1 ]o 2 ) 

where Py and r^ are the local pressure and radial distances 
of a point; Po,, , p<„ , and V„ are the pressure, density, and 
velocity of the freestream,respectively,and to is the rotational 
speed (rad/sec) of the blade. The instantaneous power 
coefficient per blade is computed from instantaneous blade 
surface pressures by 


where AA;j is the elemental surface area formed by four 
neighboring grid points, r ; j is the radial distance of the 
element, D is the diameter of the propeller, and Q is the 
rotational velocity (rev/sec). The subscript 41 denotes the 
component in the tangential direction. 

The predicted total time-mean power coefficient of the 
propeller (for two blades) for test 8 (M = 0.5, J = 3.062) is 
2 percent higher than the measured value of 0.361, 
whereas thatfortest6(M=0.2,J=0.881)is5 percenthigher 
than the measured value of 0.251. Figure 4 shows the 
variation of the single-blade power coefficient with azimuth 
angle for the two test cases considered. The blade starts at 
<(>=0, and the variation of the power coefficient is shown for 
three revolutions of the blade. The expected sinusoidal 
variation of the blade loading due to angular inflow is clearly 
observed in both cases. Theamplitudeof the power coefficient 
varies ±65 percent of the mean for test 8 and ±9 percent of 
the mean for test 6. 

A Fourier transform of the third-cycle blade power 
coefficient variation gives CP= ao* aj cos cot -b, sin cot, where 
a<),aj,and b; are the Fourier coefficients and t is the time. 
The loading spectra for tests 6 and 8 are shown in figure 5. 
The first harmonic dominates the loading, which lags the 
blade motion for tests 8 and 6 by angles of 10.5° and 17.5°, 
respectively. 

The chordwise distributions of the steady blade-surface 
pressures computed for test 8 agree well with experimental 
data (ref. 3). Figure 6 shows the unsteady blade-surface 
pressure as a function of azimuth angle for the inboard radial 
station r/R= 0.64 Ion the pressure surface for test 8. Forthe 
pressure waveforms given in this paper, 0° corresponds to 
top-dead-center as in the presentation of experimental data 
(ref. 2). The measured sinusoidal response of the propeller 
operating at an angle of attack is qualitatively reproduced by 


the computations. The phases of the predicted waveforms 
agree well with data, and the predicted levels of the unsteady 
pressure are lower than the experimental levels. The data 
show that the blade has the largest response slightly down- 
stream (x/c = 0. 1) of the leading edge and has progressively 
less response towards the trailing edge. In the predic- 
tions, the largest response occurs near the leading edge 
(x/c = 0.049), and the response reduces gradually towards 
the trailing edge. The discrepancy between measured and 
predicted maximum response locations may have resulted 
because the leading-edge blade geometry of the test is not 
reproduced exactly in the computations. 

The pressure waveforms on the suction surface of the 
blade are shown in figure 7. The predicted response and 
agreement with data are similar to those on the pressure side, 
but both the predicted and measured largest responses occur 
near the leading edge (x/c = 0.049) on the suction surface. 
Note that the pressures on the suction surface are nearly 1 80° 
out of phase with those on the pressure surface. The 
magnitudes are underpredicted as on the pressure surface. 
The pressure-time histories at the outboard radial station 
(r/R = 0.906) are shown in figure 8 for the pressure surface 
and in figure 9 for the suction surface. The agreement of the 
predicted level and phase of the pressure signals at this 
outboard radial station is similar to that at the inboard 
station: the phases of the pressures are in close agreement 
with data, whereas the magnitudes are underpredicted. 

Oilflow studies of surface streamlines of the propeller 
blade operating under takeoff conditions (with 0° inflow 
angle) show the formation of a leading edge vortex 
(refs. 9 and 10). The steady surface-pressure measure- 
ments for Modane test 6 indicate the formation of a leading 
edge vortex (ref. 1). A broad hump in the suction-side, 
chordwise pressure distribution at r/R = 0.906, seems to 
indicate that the vortex has rolled up. With the present grid 
resolution, the steady pressure distribution is predicted 
reasonably well (not shown). The broad hump on the suction 
side pressure distribution is also predicted well. 

Figure 10 shows the predicted azimuthal variation of the 
chordwise pressure coefficient distribution at r/R = 0.9 for 
test 6 with an angular inflow of 3°. Apparently, the vortex 
(indicated by the broad hump near the leading edge on the 
suction-side pressure coefficient curve) stays on the blade 
for the entire revolution. The experimental data are also 
shown in the figure. For test 6, at this radial station, only 
three pressure transducers were operating on each surface 
(as for test 8). The agreement of the predicted pressure 
coefficients on the pressure surface with data is good at all 
four azimuthal locations. On the suction surface, the pres- 
sure coefficients are under-predicted. The discrepancy is 
greatest at the azimuthal location <|> = 270°. As discussed 
later, the pressure waveforms on the suction surface are 
dominated by nonlinear viscous effects which are not 
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represented by the present inviscid calculations. Also note 
that the data points in this figure were derived from the 
computed-time mean pressure at each transducer location. 
(The time-mean pressures of the experiment were not 
available.) 

As indicated by figure 10, only slight changes in the size 
and strength of the vortex occur as the azimuthal location of 
the blade is varied. These changes are illustrated in figure 11. 
The figure shows the suction surface static pressure contours 
which indicate the presence of the vortex for the upward- 
moving blade (lightly loaded) and downward-moving blade 
(heavily loaded). In the present case, the closed contours 
distinguish the vortex from a low-pressure region produced 
by the flow expansion at the leading edge due to the high 
incidence angle. The variations of the size and strength of 
the vortex with azimuthal location of the blade are clearly 
shown although the magnitudes are small. 

Numerical solutions of the Euler equations predict the 
leading edge vortices on swept wings and delta wings fairly 
well (refs. 1 1 and 12) although the mechanism responsible 
for the formation of the vortex in these calculations is not 
clear (refs. 13 and 14). However, the ability of the unsteady 
Euler solutions to predict the time dependence of the vortex 
as in the present case needs to be studied further and 
compared with data. Such detailed data defining the time 
dependence of the vortex are not currently available. 

The pressure waveforms for test 6 with 3° angular inflow 
are shown in figure 12 for the pressure side at the inboard 
radial station (r/R = 0.64 1). The blade-pressure response is 
very low, and the predictions agree with the data. On the 
suction surface (fig. 1 3), the data show nonsinusoidal response 
at all transducer locations. Such a nonlinear response is not 
predicted by the Euler solution. At the outboard radial 
station (r/R = 0.906), the pressure waveforms on the pressure 
side (fig. 14) show sinusoidal response at all three chordwise 
locations, and the predictions agree well with the data. 
Figure 15 shows the suction-side waveforms for this radial 
station. Here, the experimental data show a highly nonlinear 
response across the blade chord at all three chordwise 
locations: x/c= 0.299, 0.565, and 0.698. Since thenonlinear 
response occurs over the entire chord, it may not be due only 
to the leading edge vortex. The waveform's double hump 
may indicate the formation and the convection and/or break- 
down of a separation bubble (spanning the blade chord ?) 
twice during a revolution - once during the upstroke (the 
large hump) and once during the downstroke (the small 
hump). The 9-ft-diameter propeller rotates at 1671 rpm, and 
the rotational velocity at the 90-percent radius is four times 
the axial velocity. The time per revolution is an order of 
magnitude higher than the time required for a disturbance to 
traverse the chord. This may partly explain the nearly 
inphase nature of the measured waveforms at the three 
chordwise locations. It is not clear if the phenomenon is 


related to the dynamic stall of an oscillating thin airfoil 
where a small change in incidence angle produces a signif- 
icant change over most of the chord (ref. 15). The present 
analysis does not capture the complex nonlinear response, 
which seems to result from unsteady viscous flow separation 
phenomena. Perhaps, a Navier-Stokes analysis accounting 
for viscous and turbulence effects may be able to simulate 
such nonlinear pressure responses. 


Concluding Remarks 

Numerical solutions of the unsteady, three-dimensional 
Euler equations were obtained for angular flow through a 
propeller for Mach numbers 0.5 and 0.2. For an inflow 
Mach number of 0.5, the phase of the predicted pressure 
response agrees well with measurements, whereas the 
magnitudes are underpredicted. At takeoff conditions 
(Mach 0.2) with high loading (blade setting angle at 
75-percent radius is 3 1 .6°), the prediction shows the forma- 
tion of a leading edge vortex, which is in qualitative agree- 
ment with the data. However, the highly nonlinear response 
measured over the entire blade chord on the suction sur- 
face at this Mach number is not predicted by the present 
inviscid analysis. Accounting for the viscous and turbulence 
effects with appropriate resolution of the length and time 
scales may shed light on the highly nonlinear measured 
response. 
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Figure 1 . — SR7L propfan installed in the S1-MA transonic test section. 
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Outboard station, r/R = 0.91 


Normalized Suction side ^x = 49.097 in. 
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(b) Streamwise surface. 



(d) Distribution of grid points on and near the blade surface. 
Figure 3. — Typical computational grid. 
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Figure 4. — Power per blade variation with azimuth angle. 


Figure 5. — Loading spectrum 
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Figure 6— Azimuthal variation of blade surface pressure on the pressure surface at radius ratio, r/R, 0.641 . Test 8; 
Mach 0,5: advance ratio. J. 3.062. 
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Figure 7.— Azimuthal variation of blade surface pressure on the suction surface at radius ratio, r/R, 0.641. Test 8; 
Mach 0.5; advance ratio, J, 3.062. 
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Figure 8. — Azimuthal variation of blade surface 
pressure on the pressure surface at radius ratio, 
r/R, 0.906. Test 8; Mach 0.5; advance ratio, J, 
3.062. 
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Figure 9.— Azimuthal variation of blade surface 
pressure on the suction surface at radius ratio, 
r/R, 0.906. Test 8; Mach 0.5; advance ratio, J, 
3.062. 





9 



Pressure coefficient, Cp 






Fraction of chord, x/c 

Figure 10.— Predicted and measured azimuthal variation of chordwise pressure coefficient at radius 
ratio, r/R, 0.9. Test 6; Mach 0.2; advance ratio, J, 0.881 . 
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Figure 1 2.— Azimuthal variation of blade surface pressure on the pressure surface at radius ratio, r/R, 0.641 . 
Test 6; Mach 0.2; advance ratio, J, 0.881 . 
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Figure 13— Azimuthal variation of blade surface pressure on the suction surface at radius ratio, r/R, 0.641 . 
Test 6; Mach 0.2; advance ratio, J, 0.881 . 
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Figure 14. — Azimuthal variation of blade surface 
pressure on the pressure surface at radius ratio, 
r/R, 0.906. Test 6; Mach 0.2; advance ratio, J, 
0.881. 


Figure 15. — Azimuthal variation of blade surface 
pressure on the suction surface at radius ratio, 
r/R, 0.906. Test 6; Mach 0.2; advance ratio, J, 
0.881. 
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